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General linear matrix model, Minkowski spacetime and the Standard Model 
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The Hermitian matrix model with general linear symmetry is argued to decouple into a finite 
unitary matrix model that contains metastable multidimensional lattice configurations and a fermion 
determinant. The simplest metastable state groups a two component spinor from the matrix and 
locally describes a Hermitian Weyl kinetic operator of either handedness on a hypercubic 3+1 
dimensional lattice with general nonlocal interactions. The Hermiticity produces 16 effective Weyl 
fermions at long distances by species doubling, 8 left- and 8 right-handed. These are identified 
with a Standard Model generation. Only local non-anomalous gauge fields within the soup of 
general fluctuations can produce effective massless propagators and survive at long distances, and 
the degrees of freedom to support non-anomalous gauge field subgroups of U(8)lXU(8)r are present. 
Standard Model gauge symmetries associate with particular species symmetries, for example change 
of QCD color associates with permutation of doubling status amongst space directions. Vierbein 
gravity is probably also generated. Low fermion current masses can arise from chiral symmetry 
breaking solutions of the fermion self-energy Schwinger-Dyson equations, generating W mass and 
composite Higgs states, similar to a scheme proposed by Gribov. Specific higher dimensional lattices 
with larger spinors are potentially stable but produce non-Riemannian spaces without conserved 
quadratic distances. However if the extra dimensions are compactified, the Minowski space persists 
at low energy accompanied by SM generations, potentially doubled further by duplicate zero modes 
in the compact directions to generate dark matter. The model is conjectured to have an origin in 
infinite dimensional conformal invariance and the concept of Bare Particulars. 
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I. INTRODUCTION 

If there is a fundamental physical theory it is expected 
to be based on a simple and metaphysically compelling 
principle, while at the same time providing an account 
of the numerous very specific features we observe in low 
energy physical laws. By anthropic selection, some but 
hopefully not all of these features could be accounted 
for contingently as those of a habitable stable configura- 
tion in a landscape of other stable configurations, some 
habitable and some barren. According to string theories 
that landscape is rich and highly varied, with our 3+1 
dimensional configuration of interacting Standard Model 
fermion generations not obviously typical or naturally 
selected from amongst almost countless alternative vari- 
eties. 

Such a scenario is not strictly objectionable scientifi- 
cally but there has been criticism on the grounds that too 
much anthropic freedom would render the fundamental 
theory unfalsifiable. There may be an alternative model 
to string theories that is simpler and relies less on an- 
thropic selection, more directly pointing to our universe 
amongst a set of not very dissimilar alternative configu- 
rations. The purpose of this paper is to argue that the 
simplest most symmetric of all Hermitian matrix models, 
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which may have an origin in infinite dimensional con- 
formal invariance and the metaphysical concept of Bare 
Particulars, could provide this alternative. 

Hermitian matrices that represent dynamic connection 
between all points in a discrete space are a promising de- 
vice to model the emergence of locally connected spaces. 
While not providing the sought for unification in string 
theory, they have been used to model the dynamical 
emergence of a world sheet in D < 2, and in a number of 
ways to model target space dimensions 1_. Independently 
of string theories, matrix models have been recently used 
to study spontaneous phase transitions into local multi- 
dimensional spaces [21 E]. 

In this paper the matrix rather than the world sheet 
is followed as the key element of interest in building a 
fundamental theory and, in keeping with the expectation 
of simplicity, the object of study is the unique Hermitian 
matrix model having general linear symmetry. Its scale 
invariant pathology decouples leaving a well-defined sub- 
matrix system with unitary symmetry operating on half 
the number of points, dominated by first order difference 
operators. When the total number of points N is odd 
the submatrix system has a determinant that represents 
strongly coupled nonlocal fermions. Statistical weight ar- 
guments suggest that only certain multidimensional con- 
figurations potentially are metastable spacetimes with 
quadratic distance measures. The dimensionality, sig- 
nature, particle spectrum, and effective interactions are 
tightly constrained and all look very like our Universe. 



II. THE GENERAL LINEAR MATRIX MODEL 

A. Construction 

If A is an N by N general complex matrix and H is 
Hermitian, then the transformation 



H' = A^HA 



(1) 



preserves Hermiticity. Physics is proposed to be repre- 
sented by a Hermitian matrix model invariant under (1) 
for invertiblc A, with N cither infinite or large enough to 
provide all the points needed for a discrete Planck scale 
spacetime. 

The flat integration measure in terms of real and imag- 
inary components is 



DH — Y\_ dH iiR J^J dHij R dHiji. 



(2) 



l<3 



Under (1) this transforms as DH' = det(AU) A ' DH. 
This can be seen by considering that the general complex 
matrix M = H + %A where A is also Hermitian has Ja- 
cobean J 2 where J is the Jacobean for H' or A' . The Ja- 
cobean for M' is easily calculated by multiplying M first 
on the right by A and then on the left by A^ contributing 



det(A'''A) Ar at each stage. The invariant measure DH is 
therefore DH/detH N . 

The unique and seemingly barren partition function is 
formed out of just the invariant measure which we call 
H- theory 



DH 



DH 

det H N 



(3) 



Finite effective subsystem 



We take the standard path of examining the eigen- 
values of the configurations that dominate the partition 
function, and since there is a unitary symmetry only the 
eigenvalues control the weight. 

Using the techniques of random matrix theory(3HZ] H 
is diagonalized to extract the eigenvalues to obtain 



Z H = V DXs(X) N e- v ^ 



(4) 



where V = J DU is the constant volume of the unitary 
group, s(A) is the product of the signs of the eigenvalues 
and 



V(X) 



N 



N 



N 

i<j=i 



x 3 f 



(5) 



The value of the eigenvalue integral, having a scale 
invariance, is not well defined, but nonetheless we can 
show that a well defined system consistently decouples. 

The commonest or highest weight eigenvalue configu- 
rations can be found by minimizing V . V is analogous 
to the potential energy (no kinetic energy) in a classi- 
cal one dimensional Dyson gas of N mutually repelling 
point particles of position coordinates A„ and charge +1 
attracted by a fixed nucleus at A = of charge —N/2. 
The mutual repulsion coming from the second term in V 
comes from the measure DH due to the Van der Monde 
determinant Jacobean. This is universal in Hermitian 
matrix models and normally leads to constant low-lying 
eigenvalue density for all V that are well-behaved at zero, 
unlike this case. 

The singular behaviour at zero can be shown to decou- 
ple by using the gas analogy. There is a balance between 
the attractive power of the pole and the repulsive power 
of the random statistics, so that only half the eigenvalues 
collapse. If N is even, N/2 eigenvalues will collapse into 
the nucleus to form an infinitessimal neutral atom and if 
N is odd (N+l)/2 eigenvalues will collapse to form an in- 
finitessimal ion of charge + 4 . The remaining eigenvalues 
will interact among one another as a free gas either with- 
out effect from the neutralized atom (N even) or under 
the milder repulsive influence of the central +1/2 ionic 
charge (N odd). 
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In either case the particular infinitessimal values of the 
bound eigenvalues which collapse to zero have no effect 
on the free eigenvalues, the only effect being the central 
charge neutralization. The free eigenvalues can only see 
the aggregate central charge. The divergent dynamics of 
the central atom therefore decouples and the free eigen- 
values are part of an independent submatrix system. If 
N is odd, where there is a residual central repulsive ion, 
the remaining N' = (N — l)/2 eigenvalues correspond to 
the system J DH det H, the residual determinant being 
suggestive of a pregeometric system of strongly coupled 
fermions. If N is even, where the central atom is neutral 
and invisible, the remaining N' = N/2 eigenvalues corre- 
spond to the system J DH which unlike the former case 
is empty of fermionic matter. 

Before proceeding we need to go back and make a scale 
transformation in the eigenvalue integral. The system 
j DH det H is dominated by the free gas of eigenvalues 
that now fly off to infinity under unconstrained mutual 
repulsion. However this infinity can also be decoupled by 
scaling it down into the central nucleus infinity by noting 
that the eigenvalue integral is of the special form 



lim / DXF(X) 



(6) 



where DXF(X) is scale invariant. Rescaling Aj — > Xis/p 
with s finite leads to 



this term has negligible effect on the overall shape of the 
eigenvalue density, which is swamped except for these 
near-zero eigenvalues (representing cosmic distances if 
distance is inverse eigenvalue) by the mutual repulsions 
amongst the gas of N' charge +1 particles. Consequently 
we can neglect this term in determining the dominant 
eigenvalue density, which is therefore controlled solely by 
the measure, i.e. the statistics of random configurations, 
and the confining box. 

It is argued later that the determinant creates the 
fermions of the Standard Model and their strong non- 
local interaction condenses to local propagating gauge 
and gravitational fields. The negligible local effect of this 
determinant at the pregeometrical stage shows that the 
statistical influence of the matrix integral that creates the 
local spacetime background in this model is a far stronger 
force than the forces we know about. 

If the eigenvalues are labelled X(x) in the continuum 
limit with x a real number ranging lowest to highest from 
to 1, the dominant density of eigenvalues p(X) — dx/dX 
can be determined. The force —dV/dXi on each eigen- 
value i due to the others must balance, and so the eigen- 
value density near the box boundary must go to infinity 
due to the force from one direction. We guess the simple 
form p(X) oc (1 — A 2 ) -1 / 2 in analogy with the Wigner 
semicircle law (a 2 — A 2 ) 1 / 2 for the Gaussian confining 
potential. This is confirmed by writing the continuum 
version of the force balance equation as 



lim / DXF(X) 



DXF(X) 



(7) 



Essentially the integration limit is arbitrary due to the 
scale invariance. Choosing s—1 the previous analysis is 
repeated. The well-behaved subsystem for odd N' has 
Dyson gas potential within the box — 1 < \ < 1 of 



P 



dx TTV 



= 



(10) 



Vy G (0, 1) where P is the Cauchy Principal value, then 
showing by contour integration that A = — cos(7rx) is the 
solution. See Figure [T] 



N' N' 

V ^ = - 2 E lo § A ? - E lo §( A * - x tf ( g ) 

i— 1 i<j — l 

which, since the function X 2B converges to the box 
function for large integer B, is generated by the reduced 
matrix model as B — > oo 

Z R = J DH det He~ trH2B . (9) 

The regulator reduces the symmetry from unimodular 
to unitary. 

C. Dominant Eigenvalue density 

The residual +1/2 central charge term in V from the 
determinant prevents zero eigenvalues and affects the 
first few eigenvalues neighbouring zero, but at large N' 



D. Local Difference Operator Form 

The constant eigenvalue density near zero is a uni- 
versal feature of Hermitian matrix measures and causes 
the dominant configurations, in the absence of special 
or tuned potentials, to be unitarily equivalent to a local 
first order difference operator. If H is expanded about its 
background value as H = Hb + Hi then by the unitary 
symmetry H — > U HU^ , 

det(J3> + H B ) = det (JZj + UH B U^) (11) 

and so all unitarily equivalent Hb describe the same 
physics. The local difference operator will be always 
be the one from which the physical law is written, so 
we perform the unitary transformation of re-ordering 
the eigenvalues away from their lowest to highest or- 
der to form a momentum space operator that has the 
period of a full sinusoidal cycle as shown in Figure [2] 
which for smoothly spaced eigenvalues corresponds to 
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FIG. 1. Dominant eigenvalues for the box potential. A: 
eigenvalue density ; B: eigenvalues ordered lowest to highest 



the Hermitian nearest-neighbour first order derivative 
iA. xy = | (f^x+i — $x,y+i)- Interestingly in this case there 
are no further modes because the shape of the momentum 
space eigenvalue function is precisely sinusoidal. This is 
different from other matrix models, which like the Gaus- 
sian case have decreasing components of third-, fifth- and 
higher odd nearest-neighbour terms, due to the difference 
of their eigenvalue density from the form that gives pure 
sinusoidal shape. 

Throughout we ignore overall normalisation constants 
in Hb since the constraint on the eigenvalues of Hi + Hb 
is unambiguous. 

The configuration that minimises the Van der Monde 
term will have smoothly spaced eigenvalues and be pre- 
cisely the first order difference operator on a one dimen- 
sional nearest-neighbour lattice. 

Typical configurations of H when N is large (i.e. com- 
mon configurations chosen at random from H) will still 
be approximately the nearest neighbour first order differ- 
ence operator but with random weight between adjacent 
points that reflects a random spacing s between adjacent 
eigenvalues roughly according to the Wigner surmise [7 



FIG. 2. Unitary transformation to local difference operator 
- A: lowest to highest order showing every second eigenvalue 
re-ordered B: sinusoidal result 



probability distribution se s . 

Neither the minimizing nor typical configurations of H 
are a multidimensional spacetime, but this is irrelevant to 
whether physics is described by the general linear model. 
We find ourselves in a highly untypical place on the sur- 
face of Earth because the vastly more typical places in 
outer space are inhospitable, and by a weak version of the 
anthropic principle the same can apply to the configura- 
tion of a habitable potentially stable multidimensional 
space-time in H which is otherwise dominated by the 
uninhabitable one-dimensional class. Since the matrix 
connects all points arbitrarily, any Hermitian multidi- 
mensional matrices are configurations within the system. 
Our Universe with its low energy density after the Big 
Bang can be represented by a decoherent region around 
a multidimensional configuration as long as the configu- 
ration is stable against small local fluctuations, and the 
local fluctuations are all that survive at low energy. 
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III. SPINOR SPACETIME LATTICES 

A. Metastability Criteria 

A definitive demonstration of metastability of any 
particular multidimensional configuration is not within 
grasp, but criteria based on statistical weight can be 
given. 

Since the first order difference operator lattice is stabi- 
lized by the mutual repulsion of eigenvalues in the mea- 
sure DH for a large Hermitian matrix, the same statis- 
tical forces apply to any large Hermitian submatrix of 
H. Consequently in any metastable configuration, every 
point in the vector space will have local neighbour con- 
nections as part of one or more locally stable first order 
difference operators. The larger the identifiable subma- 
trix, the more stable and the greater statistical penalty 
to be paid by disrupting the point from its neighbour 
connnections. 

Consequently a spacetime with simple regular connect- 
edness such as a hypercubic lattice threaded at every 
point with first order difference operators inherited from 
the one dimensional configuration is favoured over ran- 
dom connectedness that will strand many points from 
participation in any local difference operator. For a D di- 
mensional spacetime, each point on the lattice is passed 
through in D directions by D individually stable first 
order difference operators and will participate indepen- 
dently in crossing local submatrices. For this to occur, 
the D operators must be independent at each spacetime 
point which requires that a number of points from the 
vector space are gathered together to form an internal 
spinor degree of freedom. 

While a submatrix might not necessarily have its eigen- 
values bound by the box potential, there will be an ef- 
fective confinement provided by the environment of the 
matrix and its overall box potential bound. Consequently 
there may be departure from the pure nearest neighbour 
difference operator of the one-dimensional configurations 
and inclusion of third, fifth and greater neighbours, but 
nonetheless a local difference operator is expected. We 
also expect an element of variation in the strength of each 
link, in particular since the eigenvalues of an exactly reg- 
ular multidimensional lattice are have degenerate points 
associated with permutation symmetries, forbidden by 
the eigenvalue repulsion. 

One major consequence of the submatrix argument is 
that large zero Hermitian submatrices are highly unsta- 
ble. This exclusion criterion powerfully restricts the po- 
tentially stable configurations. 

B. 4-dimensional Minkowski Weyl lattice 

In what follows spinor spacetimes will be considered 
with total numbers of configuration space points N' even, 
which is necessary since spinors have two points per 
spacetime point. A single odd spacetime point, which 



FIG. 3. submatrices that contribute to stability 

must be part of the total N' system if N' is odd, can 
always be added at a boundary. 

We begin with the simple Hermitian hypercubic lattice 
representation of the left-handed Minkowski Weyl kinetic 
operator in four dimensions er^iA^, which is clearly a 
configuration in H since any H can be decomposed into 
four Hermitian matrices by H = d^H^. 

Assuming each dimension has L points, each matrix 
A M contains L 3 distinct one dimensional difference oper- 
ators parallel to the /i axis, one for each for each coordi- 
nate value of the other three dimensions, e.g 

A-0 X y = S Xl y 1 S X2 y 2 S X3 y 3 A Xg y g (12) 

There is no spacetime point that is not connected into 
such a difference operator. Any removal of a neighbour 
link will incur a large statistical penalty in one of the 
large stable one dimensional difference operators extend- 
ing along one of the axes. As shown in Fig. [3] the sub- 
matrices that influence the stability of the lattice are not 
just straight sets. There are submatrices linking a set 
of points that turn one or more corners, where it can be 
shown (at least for the nearest neighbour case) by trans- 
formations on the spinor matrices that the eigenvalues 
associated with such crooked sets are the same as the 
straight sets. These crooked sets help to favour the si- 
multaneous diagonalisation of all the that makes all 
dimensions simultaneosly local (e.g. [<9^,<9„] = 0) and 
also to favour the same link length in each dimension. 

From this example it follows that two or three space- 
time dimensions is highly unstable, since we still need a 
spinor index and the spinor matrix will be unsaturated, 
with at least one of the or a linear combination equal 
to zero. This does not however rule out toroidal compact- 
ification of the four dimensional case. Similarly other un- 
saturated spinor matrices such as cr M © cr M necessary for 
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a fundamental Dirac particle or any other unsaturated 
multiplicity are also unstable. 

Euclidean space Weyl operators are also not allowed 
because they are non-Hermitian. 

In summary, the Minkowski right or left-handed 3+1 
dimensional Weyl class is the only potentially stable 
fermion operator in four dimensions, and three and two 
dimensions are unstable. 

Nielsen and Rugh[51 |U] have previously made related 
general observations about the special dynamical stabil- 
ity of the 0+1 and 3+1 dimensional Weyl operator. 

The way in which multidimensional spacetime is 
treated here, as a metastable lattice configuration of 
threaded first order difference operators stabilized by the 
quantum measure DH, contrasts with other various ma- 
trix models previously applied to string theory that have 
been specifically designed with a "time" coordinate M(t) 
and a quantum mechanical action to generate an extra di- 
mension from the matrix eigenvalues [101 ITT) or that gen- 
erate a fuzzy sphere using a partition function over multi- 
ple matrices weighted by an SO(3) invariant action[2l[3]. 

The key assertion is that there is enough localising 
power in the measure alone to build a metastable mul- 
tidimensional lattice, which is the simplest potentially 
metastable configuration to potentially contain life and 
happens to be Minkowski space. The analysis that fol- 
lows argues that this configuration might actually contain 
the Standard Model and gravity 



C. 16 Fermions from Species Doubling 

For a general background configuration Hb the deter- 
minant can be converted to an action with Grassmann 
spinor fields ip ax and ip ax and the reduced system ex- 
panded about Hb up to a multiplicative constant is 

Z R = J DH i J Di/jDi/je^ {HB+Hl)i ' (13) 

where the B label on the integral records the regulator 
that confines the eigenvalues of Hb + Hj , and i is in the 
exponent to reflect the conventional format, discarding 
the consequent i~ N multiplier out the front of Z. Do- 
ing a perturbation expansion about Hi — when Hb is 
the Minkowski lattice Weyl operator produces the lattice 
Weyl equation ia^A^ip = from the stationary point. 

Because the Minkowski Weyl operator on the lattice is 
local, first order and Hermitian[rj3], it produces species 
doubling so that at the stationary point in the long dis- 
tance regime there are actually 8 left handed and 8 right 
handed species. In momentum space it is the zero values 
of sinj?^ at p^ = irk^ = ±ir as well as at p^ = that 
are necessary to provide the local derivative (see Fig [2j3 
above) that result in 2 4 states. Doubling has previously 
been seen as a problem for lattice simulations that add 
unwanted degrees of freedom, but in this specific case, 



where there are no bare gauge propagators, the extra de- 
grees of freedom can be the various quarks and leptons 
of a generation. 

Species doubling is not affected by adding 3rd, 5th 
or higher order modes in each dimension, as they just 
change the shape away from pure sinusoidal while leaving 
the zeroes unchanged, as in the Gaussian matrix model, 
nor is it affected by keeping the hypercubic connectedness 
but adding randomness onto the link weights [T3]. It has 
also been noted in the fcrmionic sector of string matrix 
mo dels [H]. 

The extra species can be absent in lattices where the 
points arc randomly connected |13l[T5] . and there are also 
specially designed variations to the hypercubic lattice 
derivative that can eliminate the doubling speciesfTS]. 
useful in computations in the lattice approximation to 
QCD. Neither of these apply to this necessarily spinor- 
dense hypercubic lattice. 

The species are best expressed for the current purposes 
in position space. We write the Hj = action as 

So = '^2-ip x ia> Jl ^(8 x .y + ( l - S ytX+ p,)^ y . (14) 

xy 

where fl is a unit lattice displacement in the +/x direc- 
tion. 

The extra species in position space are rapidly vary- 
ing in one or more direction - e.g. (— l) V2 Xy solves the 
stationarity criterion if Xy is a solution of 

_ S y , x+fi ) X y = (15) 

V 

which means that x' — a X obeys the right handed 
Weyl equation ia^A^x' — 0- Consequently the low en- 
ergy/momentum region contains not just slowly varying 
ij)y but also rapidly varying regions such as (— l) V2 Xy, 
where Xy is slowly varying. 

We can express the fermion field in terms of slowly 
varying species components ip x and l tp x where d = 
(da,di,d 2l d 3 ) by 

^ = £(-1)^/^ i ^ = E^(- 1 ) d " ( 16 ) 

d d 

where each d r takes the value or 1, d ■ x = ^ r d r x r 
and fd and fd are spinor matrix factors to be deter- 
mined. The slowly varying components need to be de- 
fined in terms of general functions of x convoluted with a 
smoothing function that confines the momenta to the re- 
gion [— 7r/2, tt/2], but there is no need for this detail here. 
Referring to the action it is easily seen using identities 
such as 

(-i^^C-i)^ = (-i) {d+d '> x (-i)^ d '' s ^s x , y+fl 

(17) 
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that the oscillatory component vanishes only when d = 
dl and for slowly varying functions the cross terms linking 
xp d with ip d for d ^ d' are therefore cancelled to zero and 
the continuum Hj = action is 

S a = J d 4 xJ2i> d fdi<t fl (-l) d ' S »- s d lJt M d . (18) 

d 

Choosing f d = i d ° (a 1 )^ (<r 2 ) d2 (er 3 )^ and f d = 
(cT 3 ) d3 (<7 2 ) <:!2 ((7 1 ) dl i do casts the spinor matrices into the 
standard form. For example, in the above example 
X = ip 0010 and fd = fd = o -2 - Note that in general 
fd fl due to the i da terms. The final canonical contin- 
uum form is 



5 = J d i xY J ^ dR i^d^ dR + ^^1^3^ (19) 

d R d L 

where the dL have an even number of ones and the dn 
an odd number. 



D. Interspecies Interactions 

Interactions between doubling species have been rarely 
considered before and only in the context of a lattice de- 
fined with a plaquette action, i.e. a bare gauge kinetic 
term, in which case the species interactions are off-shell 
artifacts. In this case, there are no bare gauge kinetic 
terms and the interactions actually create the mass shell 
of the effective gauge propagators. Rapidly varying Hi 
fluctuations produce couplings between the species at low 
energy if the rapid variation in Hi cancels the rapid vari- 
ations from the cross-terms in the species, which will be 
shown explicitly in the next subsection. 

In general there are a huge number of interactions in 
Hi, both local and nonlocal, most of which are neither 
Lorentz nor gauge invariant. A fluctuation can only be- 
come an observable low energy held if there is a long dis- 
tance correlation generated in the effective action i.e. a 
kinetic term generated from fcrmion loops that has a low 
or zero mass. Such terms can only arise for fluctuations 
possessing a gauge invariance that through Ward Iden- 
tites protects against masses of order of the lattice scale. 
It is certainly true that if we start with two bare fields, 
one sharing a gauge invariance with the fermion field and 
one not, the gauge invariant field is the only one that can 
have low mass and be seen at low energy. We extrapo- 
late from that clearcut case to assert that the soup that 
is Hi will select out at low energy a "fine tuned" subset 
of linear combinations of the Hi as an anomaly free set 
of propagating gauge fields. 

The first consequence is to ensure the metastability 
criterion requiring that only local fluctuations propagate 
- since the fermions are local, so will be the gauge fields 
since only local gauge invariances can be supported with 



a local fermion and all the vast number of non-local fluc- 
tuations are not observed at low energy. Therefore the ef- 
fective low energy interaction among the species, if there 
is any, can only be an anomaly-free local gauge theory. 
The only possibilities are gauge invariances that respect 
the global invariances of the background configuration, 
namely Poincare invariance and U(8)l ® U(8)r. 

The general possibility of propagating gauge fields be- 
ing emergent at low energy from a fundamental random 
dynamics with many nonlocal degrees of freedom and in- 
hereting locality from local fermion loops has been pro- 
posed previously by Nielsen and Rugh [8]. 

Scalar fields linking right and left handed species are 
probably ruled out, as there is no gauge invariance that 
can prevent a lattice-scale scalar mass, the well-known 
problem for fundamental Higgs in non-supersymmetric 
models. Consequently any non-zero light fermion and 
gauge field masses need to be generated by some form of 
dynamical spontaneous symmetry breaking. 

In the next section a specific dynamical mechanism 
that may operate to produce fermion and W masses will 
be discussed. First, it is shown that Hi does contain 
enough degrees of freedom to produce at least the Stan- 
dard Model interactions. It is not immediately obvious 
that this is so, despite its general random form, since the 
species are packed together non-trivially into the single 
spinor through doubling. 



E. Local Bare U(8) L ® U(8)r symmetry 

First it is shown that the configurations required to 
turn the global U(8)l <E> U(8)r symmetry of the nonin- 
teracting action into a bare local one are present in Hi. 
This group is anomalous and must therefore not be an 
effective theory at any scale, at least with only the 16 
fermions in the 8l and 8^ representations. We expect the 
anomalous currents in U(8)l ® U(8)r to generate masses 
of the order of the lattice scale in the corresponding in- 
duced propagators and decouple, leaving an anomaly-free 
subgroup as the effective gauge theory. We note that 
the non-anomalous Pati-Salam Model and the Standard 
Model gauge group are indeed subgroups. 

The task of establishing that these degrees of free- 
dom are present is done if Si = ipHiip contains terms 

of the form \[) dL cr ll G'f£,ij}t L and the conjugate expression 

■0x £ ct m G^ **Pt L where G dd is slowly varying, for all d, d! , 
and similarly for the right handed equivalent. 
First introduce the projection operator 

Pxv (20) 

a 

This has the property that it sends any rapidly vary- 
ing local function to zero - i.e 'Yl< x (f>x{~^) XT Pxy — 
^2 V Pxy(— l) Vr 4>y = - and it has no effect on slowly 
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varying functions. This can be used to describe fluctua- 
tions that connect two species and no others. Consider 
H dd> = A dd' + A dd'\ where 

Afy = fl(-l) d - x ^G<liP xy {~l) d '-yfl, (21) 

From the properties of P, ^jA dd ip links only 'tp d and 
i/j d , and using /'/ = /f = 1 we find (dropping the L 
label, which is assumed) 

M dd '^ = J d^x^GlUi (22) 

and so the first part of the task is done. 
For the conjugate expression, 

(A dd '% y = fd'(-l) d '- x ^G d i*S xy (-l) d -yf d (23) 

it is clear that A dd t links only ty d ' and tp . When 
= d' Q , * ne s P mor factors fd'fd' and fdf d are both ei- 
ther +1 or -1 and so multiply to 1. Accordingly, when 
do = d' all the 2 by 2 field matrices in species space 
linking ipd and ipd' with ipd and ipd/ are Hermitian i.e. 
1G M +T i G* with G, G i all real, forming the desired U(2) 
gauge invariance. However when do 7^ d , the spinor 
factors multiply to -1 and so in species space the field 
matrices are not Hermitian, with G 1 and G 2 pure imag- 
inary and the gauge invariance, although still present, is 
not unitary, being of the form Sip = ia z T l ip, Sip = —ia l T l , 
with a 1 and a 2 pure imaginary. This is a consequence 
°f / 7^ / > necessary to remove the negative sign in front 
of do in the kinetic operators for the time doublers to 
obtain the canonical form a^d^. 

Nonetheless it appears that the effective action will be 
just as if G 1 and G 2 were real. The effect of the pure 
imaginary form is to multiply the G 1 and G 2 vertices by 
i, which is compensated by a negative sign in the effec- 
tive gauge propagators that will be generated by fermion 
loops in the effective action. A simple example of this 
equivalence can be seen with the ordinary effective ac- 
tion of QED, where substituting —> iA^ yields the 
same Green's functions between the fermions, and the 
gauge invariance on the fermion field becomes ip — ¥ e a ip 
and ip — > ijje~ a , for a real. 

The same arguments apply to the right sector, and 
so the configuration J2d>d' H dd contains the degrees of 
freedom to support an effective U(8)l ® U(8)r gauge 
invariance, and therefore also contains the gauge fields 
to support the Standard Model gauge group. 



F. The Standard Model and Symmetries 

It does not follow from the above analysis that there 
must be an effective non-anomalous gauge symmetry 



at some scale, since there is a great surfeit of non- 
propagating degrees of freedom in Hj that will affect the 
upper limit of loop integrals, indeed the distillation into 
separate species itself dissolves at the upper limit, and 
the possibilty that there are no effective interactions at 
all has not been ruled out. Nonetheless the existence of 
the necessary degrees of freedom means that the Stan- 
dard Model, which is an anomaly free subgroup, is one 
candidate. 

There are some correspondences with symmetries of 
the species suggesting that the Standard Model subgroup 
is particularly suited to the left handed lattice Weyl oper- 
ator, compared with alternative anomaly free subgroups 
such as Pati-Salam. First, it seems reasonable that there 
should be a chiral asymmetry in the effective action, with 
the mass of Wr greater than Wl , considering that Hb is 
fundamentally left handed and the left-handed configu- 
rations in Hi will typically be more symmetrical. As an 
example of a highly symmetrical operator that is chirally 
selective, the projection operator P selects out only v L : 
i.e. at low energy ipPip = Vi,v L . 

Table|l]shows an assignment of quarks and leptons that 
reflect further correspondences. 



TABLE I. One possible set of quark and lepton assignments 





do 


di 


d 2 


d 3 




do 


di 


di 


d 3 


















1 











eh 


1 


1 


1 


1 


en 





1 


1 


1 


ULr 





1 


1 





UR r 


1 


1 


1 





ULb 





1 





1 


URb 


1 


1 





1 


u Lg 








1 


1 


URg 


1 





1 


1 




1 








1 


dRr 











1 




1 





1 





dRb 








1 





d L g 


1 


1 








d R g 





1 









Full reversal of doubling status in all dimensions is an 
especially simple transformation, and with the above as- 
signments reversing L(R) doubling status separately cor- 
responds to change of weak L(R) charge. For example, 
u L b (d=0101) reverses to d Lb (d=1010). 

Further, the discrete reflection and permutation sym- 
metries of the Weyl lattice are cubic, not hypercubic due 
to the difference between tr° and a % that gives the Lorentz 
50(3, 1) invariance of the stationary states. This means 
that species interactions relating to discrete space dimen- 
sion symmetries might be expected to survive preferen- 
tially. Change of SU(3) color in Table [I] can be seen 
to be equivalent to permuting doubling status amongst 
the space directions only. For example, changing red to 
blue corresponds in all quark fields to swapping the dou- 
bling status in the 2 and 3 directions. The electrons and 
neutrinos, having the same doubling status in all space 
directions, are singlets under this operation as required. 
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G. Effective Gauge Propagators and Mass 
Generation 

What has so far been shown is that the simplest po- 
tentially stable configuration of H is the Hermitian Weyl 
Minkowski space lattice, that the low energy limit has the 
same fermion spectrum as the extended Standard Model, 
that the fluctuations of Hi containenough degrees of free- 
dom for the effective gauge group to be a subgroup of 
V (8) l®U (8) Ri and that the Standard Model symmetries 
seem to fit neatly with lattice and species symmetries. 

The low energy effective theory will contain the only 
interactions that can survive at long distance by gener- 
ating massless or near-massless propagators, such as the 
anomaly free set of local gauge invariant interactions dis- 
cussed above. 

Finding the stationary point of the effective action is 
equivalent to solving the Schwinger-Dyson equations for 
the dressed propagators and vertices. The dressed quark 
and lepton propagators in general can acquire mass and 
gauge propagators can be generated due to fermion loops. 

In a lattice theory with bare gauge kinetic term j3F 2 
the full dressed gauge propagator D can be expressed in 
terms of the full pure gauge field propagator D g (which 
goes like 1//3) and the full fermion polarisation tensor II 
as 



D = D n 



D g UD g 



DgUDgUDg 



(24) 



summing to£>~ 1 = (l — HD g )D' g 1 . In the bare strong 
coupling limit /3 — 0, all the gauge propagation is gen- 
crated by the fermions and the full inverse propagator 
is 



D~ 



n 



(25) 



which is a non-closed form integral equation for D since 
IT contains D in the full vertex. At the Planck lattice 
scale, scale invariance and therefore the renormalization 
group is a bad approximation and we expect masses and 
couplings to be dependent on the high momentum cutoff 
defined by the Planck scale. At the one-loop level the 
effective coupling strength near the Planck scale cutoff A 
will have the form 



V 



No 
3tt 



log(A 2 /g 2 



(26) 



where Nq is the number of generations, although this 
could be substantially modified in the full propagator by 
the many other non-gauge invariant and nonlocal contri- 
butions in Hj. At longer distance scales where the dis- 
creteness of the lattice is not important and all the non- 
local and other non-gauge invariant field contributions 
have fallen away, the renormalization group analysis will 
become valid and the full general non-Abelian character 
including triple gluon vertices will appear in the running 
coupling from the full vertex in II. This scheme may be 




SU(2) 



SU(3) 



Mgut Mpl q 2 



FIG. 4. Emergence of RG running couplings from the 
strongly coupled theory 



summarised in Figure [4] for the case that the Standard 
model emerges from near an intermediate GUT scale, al- 
though alternatively the inverse couplings could meet at 
the strong coupling lattice cutoff M p i . 

In QED, strong bare couplings have been studied in 
many approximations to the Schwinger-Dyson equations 
and in Euclidean lattice simulations, where it has been 
found that when the bare coupling exceeds a critical 
value a fermion mass is generated, small compared to the 
UV cutoff when near the critical point [T7H2T)] . The bare 
strong coupling limit f3 = l/e§ = in QED has been 
studied in Euclidean lattices where it has been found 
that fermion mass is generated as long as the number 
of fermion flavours is less than an apparently critical 
number of about 30(21], although another approach sug- 
gests an asymptotic behaviour upper limit [2"2~ . As in ( 26 ) 
above, when the number of flavours goes up the screen- 
ing increases, and also the dynamically generated fermion 
mass will go down. 

Although the details of the fermion loop integrals near 
the lattice scale will differ from strong coupling QED for 
our model, a qualitatively similar mechanism to strong 
QED could operate in our fermionic strong coupling 
model to generate the current quark and lepton masses. 
From this chiral symmetry breaking, W masses and com- 
posite Higgs-like resonances would result. Gribov pro- 
posed such a scheme [23], where the current masses were 
assumed to be generated from the strong U(l)y dynam- 
ics at the Landau pole and he showed how the mass of 
the W-boson and the mass of effective Higgs resonances 
are a function largely of the heaviest (top) quark mass. 
Here, all interactions are strong at the Planck scale and 
so the Landau pole is shifted to the Planck scale. 



H. Gravity 

Gravity requires identification of the configurations to 
support a vierbein field and possibly an independent con- 
nection field |24j. and generation of curvature terms in 
the effective action from fermion loops as in the gauge 
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field case. Clearly the vierbein will be on the linking 
entries such as 5 X ,3/+/} that define the flat background 
space, where the fluctuations of Hj provide the degrees 
of freedom to support an effective vierbein field around 
the background value 6%. The connection field may 
emerge in a more complex combination. 

If the observed low-energy local invariance is re- 
stricted to volume-preserving coordinate transforma- 
tions, as might arise from the unimodular nature of 
the Lorentz invariance of Hg and/or the fundamental 
length 25J provided by the lattice, unimodular gravity 
would result, which does not respond to cosmological 
constant terms and has been proposed as solution of the 
first cosmological constant problem [33H2H]- 

The manner in which the graviton propagator is pro- 
posed to emerge from the bare strong coupling theory 
of fermions is analogous to Sakharov's idea of induced 
gravity from one loop dominance, which has recently 
been shown to produce the right gravitational coupling 
strength for a fundamental Planck length cutoff when the 
fermion loops of the 3 generation Standard Model are as- 
sumed to dominate [30]. The approach taken here could 
provide the theoretical justification for such a calculation. 



IV. GENERATIONS FROM HIGHER 
DIMENSIONS 

A. Non-quadratic distances 

Higher dimensional lattices can be the source of ad- 
ditional generations. The saturation stability require- 
ment leads to higher dimensions having non quadratic 
distances. For example, increasing the number of internal 
components from two to three leads to a saturated ma- 
trix of the form X^^A^y where X° a p — 5 a p and the other 
A* are the eight SU(3) generators. Whereas the fermions 
in Minkowski space move on light cones PpPv'^ 1 ' = 
corresponding to the zero modes of H, the zero modes of 
this nine dimensional space given by det(X ti p fi ) = obey 
a cubic polynomial equation. 



B. Compactified Six component spinor lattice 

When the number of internal components is even, in- 
teresting possibilities arise with compactification. Con- 
sider for example the saturated lattice with 6 components 
and consequently 36 dimensions, which we write in terms 
of a direct product of 2 by 2 and 3 by 3 spinor matrices 
as follows: 



. yfjj — iA^xySijCr^p + iA fl A,xy^tj a a/3 ( 28 ) 



(27) 



Note that time will always be identified with the direc- 
tion corresponding to the identity matrix in spinor space. 

If all the other 32 directions have relatively few 
(or very few) points and are closed then they would have 
a radius R smaller than measured by current energies. 
Since the submatrices in the compact directions are rel- 
atively small (maybe as small as 3 by 3), it is not the 
eigenvalue repulsion for the small matrices in isolation 
that will figure in the stabilising forces but larger crooked 
submatrices previously discussed, between sets of points 
travelling down a non-compact direction, taking a turn 
into the compact direction for a few steps and then turn- 
ing back into a non-compact direction. Such submatrices 
will tend to stabilise the points in the compact direc- 
tions to have same nearest-neighbour links as those of 
the larger directions, and so the stable operator A x A^, y An 
along a compact axis will be simply the closed difference 
operator. 

The momenta (eigenvalues of the A^A xy ) will then 
be spaced in units of order 1 / R and at low energies there 
will not be enough energy to excite transitions between 
momenta in the compact directions, meaning the parti- 
tion function will be four dimensional with 3 generations 
from the first term in (281 above. 



Consider that our four dimensions are the four a = 
directions, and label the other 4 times 8 directions with 
/iA, where A ranges from 1 to 8, so that 



Recalling that the Standard Model species arise due 
to there being two zero modes in each of the large di- 
mensions, there is the question of whether the compact 
directions produce any more doubling. Each compact 
direction that has doubling will double the number of 
generations seen in Minkowski space up to a maximum 
of 2 32 = 4 billion times! For an odd number of points 
the closed difference operator has one zero eigenvalue, 
and therefore no doubling. For an even number greater 
than two there are two zero eigenvalues and therefore 
doubling. However the eigenvalues generally pair for the 
even sizes and this may destabilize them compared to 
odd sizes. 

Clearly lattices with 4,8,10,12 or higher-spinors can 
also be considered. 

The additional generations from compact direction 
doubling can have different effective interactions from 
the three a — generations, depending on the form of 
Hi in the compact direction that survives in the low en- 
ergy limit. A form without dependence on the compact 
direction will lead to the full set of interactions, and a 
form that has operators like P in the compact directions 
will decouple or alter the interactions of some or all of the 
species. Such behaviour might be needed to produce a vi- 
able dark matter candidate, since regular standard model 
Dirac neutrinos appear to be disfavoured by the data[3T] , 
and could also generate brane states with only gravity 
in the extra dimensions, or mirror dark matter [32]. We 
also note that additional anomaly-compensating species 
could allow the high energy gauge group to be expanded 
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to SU(8)l <& SU(8)r which been reported as a candi- 
date GUT that can account for the low energy coupling 
strengths without resorting to supersymmetry[3~3~j . 

Derivation of the correct shape of the propagating in- 
teractions in extra dimensions from the random dynamics 
is needed to decide amongst the possibiities. 



V. POSSIBLE UNDERLYING PRINCIPLES 
A. Infinite Dimensional Conformal Invariance 

If we accept that the observed low energy physics could 
emerge from the general linear matrix model, why would 
the fundmental theory be this matrix model? Matrix 
models rip apart the rules of a local dimensioned space- 
time manifold and allow each of the N points to dynami- 
cally link with all others. There is clearly an analogy be- 
tween the general linear transformation and a coordinate 
transformation in an N (assumed infinite) dimensional 
complex space, and between the Hermitian matrix and a 
metric or vielbein over the complex space. 

In Einstein gravity the ordered continuous local space- 
time is parametrized by four coordinates x^ and the gen- 
eral coordinate invariance x —> x'(x) confines the space- 
time dynamics to the metric g^ v {x) which is needed 
to have an invariant spacetime displacement ds 2 — 
9 l iu{x)dx tl dx v . 

Consider the same process applied to a pregeometrical 
theory of N discrete points where each point a is assigned 
a separate complex number z a and demand invariance 
under general analytical transformations z a —> z' a (z). 
This yields the transformation of the displacement dz' a = 
{dz' a / dz h )dz h = A -1 (z)?dz and its conjugate, an invari- 
ant "displacement" dz T Hdz in the N dimensional com- 
plex space is provided by a (non-positive definite) metric 
Hermitian matrix H(z, z) transforming as H — > H A, 
and and invariant 'volume' is det HDzDz. 

Generalizing the simple invariant matrix integral 
j DH of H-theory to the simple conformally invariant 
functional integral J VH over the metric matrix H (z, z) 
doesn't add any new physics, as the z and z dependence 
is a dummy label that can be ignored in the absence of 
a non-trivial integrand, taking us back to the simple ma- 
trix model. Nonetheless it does provide a "derivation" of 
H and its transformation properties as those of a metric 
required to deliver the conformal invariance. There may 
also be more complicated invariant integrals involving in- 
tegrands with actions and curvature terms that should be 
investigated. 

It is worth noting that string theory manages to de- 
duce local gauge particles and gravitons by demanding 
conformal invariance over a world sheet embedded in lo- 
cal spacetime parametrized by a single complex number 
z, and in a sense one could interpret this proposal as a 
generalization of string theory, but the interpretation is 
quite different on its face and the integrands are trivial. 



B. Bare Particulars 

The arguments of the previous section embed a funda- 
mental 'gauge' principle into a pregeometrical quantum 
space spanned by discrete entities, half of which end up 
as points of a lattice spinor spacetime, and the propagat- 
ing gauge and gravitational fields in spacetime are the 
remnants of this principle localised by matrix statistics. 

The gauge principle is more like a deductive process - 
we apply numbers to span a configuration space and then 
deduce the "cleaning" forms such as invariant integrals 
over gauge fields or metrics or matrices because we want 
the numbers (or parts of them) not to matter. Why then 
is this process of application and cleaning fundamental 
to physics? 

As an attempt at an answer, imagine that the entities 
are no more than a discrete plurality of independent en- 
tities that are maximally simple in having no properties 
at all. These are not like featureless points represented 
in a space, that actually have distinguishing spatial coor- 
dinates and relations, or like featureless points on a line, 
that actually have sequential order and all the arithmetic 
that follows, but bare entities that have no intrinsic or 
relational properties. 

Bare entities suffer from a representational complexity. 
The complexity comes about because conceptual spaces 
that we are accustomed to think of things as being 'in' 
introduce something from the outside - the coordinate - 
which is laden with properties. 

However representation must involve associating some- 
thing with the entity and expressing a form with the 
something. If we just decided that since bare things are 
featureless we won't associate anything with them, we 
would have represented nothing, not a set of featureless 
entities. Associating a set of numbers {z} to the set of 
entities inherently introduces properties, the algebraic re- 
lations (e.g. 1+2=3 or i 2 = —1), so static numbers do 
not represent bareness. Nonetheless the representational 
tool has to be algorithmically generated, i.e. some sort of 
number, and the only way of writing down a form with 
numbers that has no algebraic relations is to sum over 
all number values in such a way so that at each point 
in z a -space, z a doesn't matter. Hence the need for also 
associating the metric matrix with pairs ab having its 
transformation law, and for an invariant sum J DH over 
the metric that allows the transformation to be realised 
as a change of variable inside the sum. 

Adopting the continuous complex numbers for z a and 
conformal invariance for the symmetry makes some sense 
as the representational choice because, unlike the integers 
or reals, the complex numbers are algebraically complete, 
able to represent all (commutative, associative and dis- 
tributive) algebraic relations, written as conformal equa- 
tions F(z) = 0. One may question whether the complex 
numbers contain enough algebraic relations to apply and 
then clean off - what about noncommutative and nonas- 
sociative algebraic relations of the quaternions and octo- 
nions, or the sedenions? And should we stop at cleaning 
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the entities a represented by z a - what about thinking 
of the entity pairs ab as another composite entity and 
cleaning the matrix H a t, with another matrix -ff a f, ;c d? 

There are many uncertainties but the key suggestion 
to provide a path further forward is that the fundamental 
entities are very bare, maybe utterly bare, and they defy 
the arbitrary properties of numbers. 

There is a similarity or perhaps identity of the bare 
entity concept with the Bare Particulars long considered 
by some metaphysical philosophers to be required to rep- 
resent the fundamentally 'thingy' essence of entities that 
make them something other than abstract. The Bare 
Particular associated with an entity is the entity stripped 
bare of all the properties it instantiates [31] • Sider[3"5] has 
drawn an analogy between entities associated with rela- 
tional but not intrinsic properties and points of classical 
spacetime, and a further analogy between what he calls 
utterly Bare Particulars (entities associated with neither 
relational nor intrinsic properties) and the singletons of 
Lewis's set theory. So perhaps there might be a connec- 
tion betwen physics and the foundations of set theory, 
which itself is a foundation of arithmetic. We have al- 
ready begun to suspect that the foundations of physics 
and the foundations of mathematics are the samc[36j. 

VI. CONCLUSION 

The fundamental physical theory ought to provide an 
explanation for everything, including Minkowski space- 
time, spin, fermions with Pauli statistics, gauge in- 
variance, Standard Model generations, chirality, gravity, 
mass and even quantum field theory itself. Some of these 
features can be anthropically selected but the theory has 
to tell us which ones. 

H-theory potentially provides the explanations of these 
features in sometimes unexpected ways and arguably 
points much more strongly to our familiar laws than 
string theory, although by conformal invariancc it might 
be related. It might seem unsatisfactory that such a nu- 
merically pathological expression as J DH could be pro- 
posed as a fundamental expression, but pathology may 
be an inevitable price to pay to model bare entities using 
numbers. Alternatively there may be a better formu- 
lation that also generates a finite matrix model of the 
general form J DH det H on which all the predictions 
depend. 

In the scenario painted here, physics is modelled by 
the general linear matrix model because the fundamen- 
tal entities are Bare Particulars, requiring representation 
with complex numbers in an expression having infinite 
dimensional conformal invariance. Spacetime is not a 
consequence of the quantization of curvature that leads 



to disordered spacetime foams, but of the enormously 
strong statistics of the nonlocal Hcrmitian matrices that 
lead to ordered robust local lattices of first order differ- 
ence operators. The Minkowski Weyl 3+1 dimensional 
spinor lattice has a central role as the only potentially 
stable multidimensional configuration with quadratic dis- 
tances. Time is the difference operator associated with 1 
in spinor space. Fermionic Grassmann matter fields are 
not a consequence of supersymmetry but are a remnant 
determinant from the underlying pregeometrical symme- 
try. The 16 Weyl spinors of the Standard model gen- 
erations are unavoidable and intimately linked with the 
lattice spacetime via species doubling. Generations and 
possibly dark matter are from compactificd higher dimen- 
sions. Local gauge and gravitational fields are surviving 
low energy anomaly-free effective gauge theories from a 
soup of general random fluctuations strongly interacting 
with the local fermions. 

There are many directions for future work to explore 
the several unproven conjectures and possibilities pre- 
sented here, including: 

• rigorous derivation of the reduced system; 

• the stability criteria for lattices; 

• the random dynamics mechanism ; 

• emergence of the Standard Model gauge group; 

• fermion and W mass patterns; 

• vierbein propagator and the cosmological constant; 

• cosmological consequences of the locally flat robust 
background; 

• effective propagators in the extra dimensions; and 

• representational consequences of the bare entity hy- 
pothesis. 

Some computer simulations might be tractable al- 
though representing the full non-local set of interactions 
will severely restrict the achievable lattice size. Clearly 
a major positive outcome would be to confirm that the 
Standard Model gauge group arises as the long distance 
effective interaction between the 16 species. 
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